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Abstract 

Background: In order to get global molecular understanding of one of the most important crop diseases 
worldwide, we investigated compatible and incompatible interactions between Phytophthora infestans and potato 
{Solanum tuberosum). We used the two most field-resistant potato clones under Swedish growing conditions, which 
have the greatest known local diversity of P. infestans populations, and a reference compatible cultivar. 

Results: Quantitative label-free proteomics of 51 apoplastic secretome samples (PXD000435) in combination with 
genome-wide transcript analysis by 42 microarrays (E-MTAB-1515) were used to capture changes in protein 
abundance and gene expression at 6, 24 and 72 hours after inoculation with P. infestans. To aid mass spectrometry 
analysis we generated cultivar-specific RNA-seq data (E-MTAB-1712), which increased peptide identifications by 
17%. Components induced only during incompatible interactions, which are candidates for hypersensitive response 
initiation, include a Kunitz-like protease inhibitor, transcription factors and an RCR3-like protein. More secreted 
proteins had lower abundance in the compatible interaction compared to the incompatible interactions. Based on 
this observation and because the well-characterized effector-target CI 4 protease follows this pattern, we suggest 40 
putative effector targets. 

Conclusions: In summary, over 1 7000 transcripts and 1000 secreted proteins changed in abundance in at least one 
time point, illustrating the dynamics of plant responses to a hemibiotroph. Half of the differentially abundant 
proteins showed a corresponding change at the transcript level. Many putative hypersensitive and effector-target 
proteins were single representatives of large gene families. 

Keywords: Potato, Desiree, Sarpo Mira, SW93-1015, Secretome, Apoplast, Resistance, Hypersensitive response, 
Phytophthora infestans 



Background farmers' income and food security in developing coun- 

The oomycete Phytophthora infestans, the cause of po- tries [1]. In addition, new regulations designed to reduce 

tato (Solanum tuberosum) late blight disease, is one of threats to the environment will limit the availability of 

the most confounding plant pathogens. Despite over a fungicides, and conditions favoring overwintering spores 

century of resistance breeding, fungicide use, and other due to the effects of global warming, are expected to in- 

control measures, it is still a major threat to sustainable crease problems with this disease [2]. 

food production worldwide. P. infestans is responsible P. infestans has a hemibiotrophic life cycle [3] and bi- 

for global annual costs of at least €5.6 billion in control phasic growth, with an initial biotrophic phase and subse- 

measures and crop losses and is especially a threat to quent necrotrophic phase of infection [3,4]. The successful 

suppression of plant defense by effector molecules leads 
to a compatible interaction between plant and pathogen, 
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effectors in the P. infestans genome have been identified 
[6,7] and increasing evidence for their role in establishing 
ETS exists. However, knowledge about their targets in the 
plant host is scarce, mainly because of the limited avail- 
ability of hemibiotrophic pathogens with susceptible inter- 
actions with the model species Arabidopsis [8]. Some of 
the P. infestans extracellular effectors are inhibitors that 
target defense-related proteins such as proteases and glu- 
canases, and processes such as programmed cell death in 
plants. For example, Kazal-like extracellular serine prote- 
ase inhibitors EPI1 and EPI10 inhibit the P69B subtilisin- 
like serine protease in tomato [9,10], while others such as 
the cystatin-like protein target cysteine-proteases [11] and 
EPIC1 and EPIC2B inhibit papain-like cysteine protease 
RCR3, in the same plant [12]. These latter two inhibitors 
also target cysteine protease C14 by direct binding in the 
extracellular compartment [13]. In addition to extracellu- 
lar inhibition of C14 by EPIC1 and EPIC2, an intracellular 
RXLR effector, AVRblb2, can prevent C14 secretion to the 
apoplast [14]. On the other hand, P. infestans has several 
xyloglucan-specific endoglucanases while the potato gen- 
ome includes ten clustered genes for xyloglucans-specific 
endoglucanase inhibitor proteins (XEGIPs) [15]. 

Apart from PAMP (pathogen associated molecular pat- 
terns) recognition, plants have evolved a second mechan- 
ism of pathogen recognition through resistance genes, via 
direct or indirect interaction with effectors. Through this 
mechanism plants can avoid ETS, leading instead to an in- 
compatible interaction between plant and pathogen, re- 
ferred to as effector- triggered immunity (ETI). This is 
often accompanied by a hypersensitive response (HR) 
[5,16]. PAMP and effector recognition leads to induction 
of biotic stress signaling involving MAP kinase cascades 
and induction of defense-related genes through phosphor- 
ylation of WRKY [17] and MYB transcription factors [18]. 
WRKY8 is involved in plant basal defense in tobacco, and 
is phosphorylated by MAPKs [19]. In addition, both jas- 
monic acid and salicylic acid are required for activation of 
PAMP-induced defense responses in potato [20]. 

Genome-wide transcript profiling has revealed many 
similarities between incompatible and compatible interac- 
tions, but higher numbers of differentially expressed genes 
were found during an incompatible Hyaloperonospora- 
Arabidopsis interaction [21]. However, genes with altered 
expression exclusively in compatible interactions have also 
been identified [21-23]. Microarrays based on cDNA 
clones [24-26], and more recently DeepSAGE transcrip- 
tome sequencing analysis, [27] have developed an under- 
standing of compatible and incompatible interactions 
between P. infestans and potato, but none of these studies 
were based on the sequenced potato genome. 

Here, we present a genome-wide expression profiling 
of compatible and incompatible potato-/ 3 , infestans in- 
teractions in combination with quantitative apoplastic 



secretome analyses. Although apoplastic profiling to find 
representative protein families has been performed in 
rice, chickpea, alfalfa and Arabidopsis unaffected by 
plant pathogens [28,29], and in response to pathogens to 
a limited extent [30,31], quantitative analysis of apoplast 
proteins during the progress of any compatible and in- 
compatible interactions has not been presented previ- 
ously. The recent release of the potato genome [32], 
combined with clone-specific sequence information de- 
termined by RNA-seq, increased sensitivity of current 
mass spectrometers, and new methods for label-free 
quantitative proteomic analysis enabled our study. These 
technological advances present new possibilities to help 
understand mechanisms in plant-pathogen interactions 
as well as to identify candidates for resistance against 
pathogens. 

The Nordic countries have high P. infestans diversity, 
with both of the P. infestans mating types required for 
sexual reproduction and a suitable environment for 
sexual recombination [33]. Since these conditions pro- 
vide an excellent environment to evaluate sustainable 
resistance against P. infestans, we have chosen the two 
most resistant clones based on several years of field tri- 
als in Sweden [34]. In addition, we monitor one com- 
patible interaction, allowing us to find putative effector 
targets, many of which are single proteins from large 
gene families. 

Results 

Overview of differentially expressed genes and changes 
in protein abundance 

In order to study stress responses in plant-oomycete in- 
teractions, we performed time series analyses of two in- 
compatible interactions (Sarpo Mira and SW93-1015) and 
one compatible interaction (Desiree) with P. infestans 
(SE-03058). This was performed at a global level using 
genome-wide microarrays to determine gene expression 
changes and short ID-gel separation followed by MS/ 
MS analysis on a LC-coupled Orbitrap mass spectrom- 
eter to determine quantitative changes in apoplastic 
protein levels. Sarpo Mira and SW93-1015 have similar 
HR initiation, but later Sarpo Mira has a more typical R 
gene-mediated expansion in HR lesions [34]. For the 
compatible interaction, Desiree was chosen because it 
is the most used potato cultivar in molecular and 
physiological studies. 

In our genome-wide array, 17451 transcripts were sig- 
nificantly differentially expressed relative to uninoculated 
control in at least one time point in one or more clones 
during the course of P. infestans challenge (Table 1). This 
demonstrates that this biotic stress condition has a pro- 
found effect on transcriptional responses in potato. 

In order to improve the peptide identification by mass 
spectrometry, we generated about 52 million clean pair- 
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Table 1 Number of differentially expressed transcripts in 
the microarray experiments in potato clones Desiree, 



Sarpo Mira and SW93-1015 and the number of proteins 
with differential abundance in the three potato clones at 
6, 24 and 72 hpi 



Differentially expressed transcripts 


Potato clone 


6 hpi 




24 hpi 




72 hpi 






Up 


Down 


Up 


Down 


Up 


Down 


Desiree 


1449 


1602 


633 


618 


6009 


6109 


Sarpo Mira 


1362 


1440 


543 


351 


331 


417 


SW93-1015 


1169 


1131 


1501 


1782 


2554 


2989 


Proteins with differential abundance in the secretome 


Potato clone 


6 hpi 




24 hpi 




72 hpi 






Up 


Down 


Up 


Down 


Up 


Down 


Desiree 


218 


76 


194 


82 


289 


225 


Sarpo Mira 


307 


74 


353 


66 


375 


127 


SW93-1015 


273 


45 


397 


41 


307 


38 



end reads from RNA-seq data for each of the three po- 
tato clones. Contigs were constructed by de novo tran- 
scriptome assembly and subsequently ab initio gene 
prediction algorithms were used to predict ORFs, which 
were then combined in a customized database for pep- 
tide identification. We found more than 3000 transcripts 
(ORFs) for each potato clone that lacked BLAST hits in 
the S. tuberosum predicted proteome (Additional file 1: 
Table SI) and an overall 17% increase in peptide identifi- 
cation was observed. From the apoplast of the three 
potato clones a total of 12591 peptides were identified. 
Peptides that appeared in at least two replicates per time 
point were kept, resulting in 5055 peptides correspond- 
ing to 1639 proteins, which were used for quantitative 
analysis. Of the 1639 identified proteins, 1075 were 
found to be differentially abundant in at least one time 
point during the infection, and of these 785 had Potato 
Genome Sequencing Consortium (PGSC) annotations. 
In total, 830 of the identified proteins could be classified 
into 433 different protein families by InterProScan ana- 
lysis [34,35]. Based on TargetP [36] analysis for secretory 
signals, 517 of these proteins were predicted to have a 
secretory signal, which is in line with other apoplastic 
secretome studies [29]. Among the most prominent 
protein families were subtilase family proteins, peroxi- 
dases, protease inhibitors, GDSL-lipases, and pectin es- 
terases (Figure la). Several proteins that are classified as 
pathogenesis-related (PR) proteins were identified both 
in infected and un-infected samples in all three potato 
clones, including thaumatin, glucanases, glucosidases, 
P69 proteins, peroxidases, serine proteases, cysteine pro- 
teases, PRp27 protein, subtilases, chitinases, osmotins 
and lipid transfer proteins. Many PR protein family 
members are thus present in the apoplast already prior 



to infection and could be putative effector targets very 
early in plant-pathogen interactions. 

At 6 hours post infection (hpi) and 24 hpi, similar total 
peptide intensity was observed in both compatible and 
incompatible interactions, suggesting similar amounts of 
protein at these time points. However, at 72 hpi we ob- 
served significantly higher total peptide intensity in the 
compatible interaction (Desiree) compared with the two 
resistant potato clones (Figure lb) although fewer pro- 
teins with higher abundance were observed in Desiree 
(Table 1). At 72 hpi in Desiree, glucanases and glucosi- 
dases comprised a major part of the total apoplast pro- 
tein intensity, whereas total peptide intensity for all 
identified proteases and peptidases increased in all three 
clones (Additional file 2: Figure SI). This indicates a 
more complex response of the protein composition in 
the incompatible interactions. For proteins that showed 
a change in abundance in the apoplast during infection 
only 392, or 50%, of the corresponding transcripts 
showed differential expression (Figure lc). 

A large number of transcripts were found to be differ- 
entially expressed at 6, 24, and 72 hpi in both compat- 
ible and incompatible interactions (Table 1). Different 
sets of transcripts were differentially expressed between 
earlier time points (6 and 24 hpi) and later time point 
(72 hpi) (Figure 2a), indicating a biphasic response to 
pathogen inoculation in both incompatible and compat- 
ible interactions. 

Similar numbers of induced transcripts were observed 
in all three clones at 6 hpi (Figure 2b). However in Sarpo 
Mira fewer differentially expressed transcripts were 
observed compared with the other two potato clones at 
later time-points (Table 1, Figure 2). Uninfected SW93- 
1015 had around 4000 transcripts that were differentially 
expressed compared to un-infected Desiree controls. 
These included SA induced defense marker PR1 gene 
family members (DMG400005109; DMG400005111), and 
JA induced defense marker [37] MYC2 (DMG400017535). 
This supports our earlier conclusion that this clone has a 
paranoia trait, a weak constitutively active defense [34]. 
Furthermore, GO terms related to plant cell death and 
defense such as programmed cell death, methyl salicylate 
esterase activity, apoptotic process, and innate immune re- 
sponse were enriched in SW93-1015 control samples 
compared to Desiree control (Additional file 3: Table S2). 
SW93-1015 generally had more differentially expressed 
transcripts than Sarpo Mira at 24 and 72 hpi (Table 1), 
most of which were specific to SW93-1015 (Figure 2b). 
However, a large number of transcripts with altered ex- 
pression in SW93-1015 were also differentially expressed 
in Desiree (Figure 2b), indicating that SW93-1015 might 
recognize more cues from P. infestans than Sarpo Mira. 

Similar to the transcriptome patterns, several proteins 
that differed in abundance overlapped between 6 and 
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Figure 1 Characterization of secreted proteins, (a) Number of genes found in identified protein families in the apoplastic secretome based on 
conserved domains and motifs using InterProScan. (b) Mean total peptide intensities for each sample across the four replicates. Error bars 
indicate standard deviation of the mean, (c) Venn diagrams of differentially expressed transcripts in the microarray and differentially abundant 
proteins in the apoplast. Only the 785 proteins with Potato Genome Sequencing Consortium (PGSC) annotations from 1075 differentially 
abundant proteins were used for comparison. 



24 hpi within clones, as seen in the volcano plots 
(Figure 3a). At 72 hpi, however, several more proteins 
differentially accumulated compared with earlier time- 
points (Figure 2a). In addition, a larger number of pro- 
teins decreased in abundance in Desiree compared with 
Sarpo Mira and SW93-1015 (Figure 3a, Table 1). A corre- 
sponding difference in down-regulated transcripts was not 
observed. Instead, a much larger number of up-regulated 
transcripts was seen in Desiree at 72 hpi compared with 
the incompatible interactions, as mentioned above. A 
large number of proteins that increased in abundance 
were common in all three potato clones (Figure 3b). 
In addition, a higher percentage of proteins with in- 
creased abundance were observed during incompatible 



interactions in the apoplastic secretome data (Figure 3b) 
compared with transcript data where larger number of dif- 
ferentially expressed transcripts were observed in the 
compatible interaction (Table 1). 

Responses during incompatible interactions 

A hallmark for plant resistance towards pathogens is the 
initiation of the HR reaction, but in crops relatively little 
is known about the genes involved. Among the differen- 
tially expressed transcripts we found 92, 112, and 68 
transcripts at 6, 24, and 72 hpi, respectively, that were 
up-regulated in both resistant potato clones but not in 
susceptible Desiree (Figure 2b). Some of these genes 
such as the MYB transcription factors, glutaredoxins, 
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Figure 2 (See legend on next page.) 
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(See figure on previous page.) 

Figure 2 Transcript analysis across the time course in Desiree, Sarpo Mira and SW93-1015. (a) Volcano plots showing differentially 
expressed transcripts and overlap between the different time points. Blue dots at 24 and 72 hpi indicate transcripts with present at previous 
time points and grey dots indicate unique transcripts for that time point, (b) Venn diagrams of down-regulated transcripts (upper part) and 
up-regulated transcripts (lower part) regulated in the microarray experiment in the three potato clones at 6, 24, and 72 hpi from left to right. 



RING zinc finger proteins and U-box proteins, have 
members that have previously been associated with re- 
sistance (Additional file 4: Table S3). There were 29 dif- 
ferent GO terms uniquely enriched for the early (6 hpi) 
incompatible interaction (Figure 4). Similarly we found 
59, 131, and 66 secreted proteins at 6, 24, and 72 hpi, re- 
spectively, with increased abundance in the two incom- 
patible interactions but not in the compatible interaction 
(Figure 3b). These include a subtilisin-like protease, 
Kunitz-like protease inhibitors, lipid transfer proteins, 
defensins and strictosidine synthase. These genes and 
proteins can be regarded as candidate involved in initi- 
ation of HR and resistance against P. infestans. 

MapMan pathway analyses reveal that a majority of 
the represented MYB transcription factors are down- 
regulated at 6 hpi in Desiree, whereas the opposite is 
true for Sarpo Mira and SW93-1015 (Additional file 2: 
Figure S2). In contrast, several WRKY transcription fac- 
tors are down-regulated in all three clones at 6 hpi. 

In order to identify testable HR initiation candidates, 
we took a conservative approach and identified genes 
that were 4-fold up-regulated only in the resistant culti- 
vars at 6 or 24 h, at the RNA (Table 2) and protein 
(Table 3) level. Only four of the 25 transcripts were dif- 
ferentially regulated in the apoplast (highlighted bold in 
Table 2). The 49 secreted proteins which accumulated 
early in resistant clones indicate that this compartment 
is active in HR initiation or other early resistance mech- 
anisms (Table 3). Notable proteins from the apoplast are 
the RCR3-like Phytophthora-inhibited protease 1 (PIP1), 
the aspartic protease nepenthesin-1 and a Kunitz-like 
protease inhibitor. These findings show that the proteo- 
lytic machinery of potato plays an important role in 
defense against pathogens. In line with these findings, 
tomato plants with mutated Rcr3, which is related to 
PIP1, have enhanced susceptibility to P. infestans [12]. 
Several hydrolytic enzymes are also identified and these 
might exhibit novel catalytic specificities evolved to 
assist in ETI. 

General defense and the susceptible interaction 

There was a clear overlap between differentially 
expressed transcripts in both the incompatible and com- 
patible interactions (Figure 2b) demonstrating that there 
are common sets of genes differentially expressed inde- 
pendently of the type of interaction. At 6 hpi there is a 
large overlap between all three clones in enriched GO 
terms (Figure 4) with many GO terms linked to primary 



metabolism, but also general signaling mechanisms such 
as protein tyrosine and serine/threonine kinase activity. 
A change in primary metabolism and brassinosteroid- 
related processes is then seen throughout the time 
course in all three clones (Additional file 3: Table S2). 
Between the three clones at 6 hpi, the MapMan pathway 
analysis shows a striking similarity for C2H2-CO-like 
and C2H2-Dof zinc finger containing transcription fac- 
tors indicating that these might be involved in general 
defense (Additional file 2: Figure SI). There are a large 
number of GO terms only enriched in the compatible 
interaction (Figure 4). In hormone-related processes, 
auxin transport and stimulus GO terms are unique to 
the compatible interaction, and this is also true for gib- 
berellin biosynthesis and stimulus. At 72 hpi a large 
number of GO terms are significantly enriched in 
Desiree. It is the only time point where functional 
groups associated with jasmonic acid and systemic ac- 
quired resistance are present, and is in line with the start 
of the necrotrophic phase of the pathogen. Also, several 
histone-related genes can be seen to be up-regulated in 
Desiree, whereas this is the case for only a few members 
in SW93-1015 and none in Sarpo Mira. 

Screen for putative effector targets in the apoplast 

Although the use of effector targets in resistance breeding 
is desirable [38], very few of such proteins have been 
identified in oomycete-plant interactions. We detected 
many secreted proteins that decreased in abundances 
during the compatible interaction. This is in line with a 
suppression of some secreted proteins by pathogen effec- 
tors or a faster degradation due to protein-complex 
formation. P. infestans effector can suppress secretion of a 
defense protein is exemplified by the cysteine protease 
C14 [14]. We found that the C14 protease decreased 
sharply at 6 hpi in the Desiree apoplast even though it 
increased at transcript level at this time point validating 
Bozkurt et al's findings in tomato (Figure 5). This inspired 
two strategies to screen for putative effector targets in our 
quantitative apoplastic secretome data. 

Firstly, we selected proteins with a minimum of 8-fold 
decreased in abundance in the susceptible clone at 6 and 
24 h after infection. This selection produced 40 candi- 
dates, many of which are single proteins from large gene 
families (Table 4). These proteins may also be putative 
targets in the incompatible interactions because at early 
time points secreted effectors are expected to be present 
in all interactions. 
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Figure 3 (See legend on next page.) 
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Figure 3 Differentially abundant apoplastic proteins across the time course in Desiree, Sarpo Mira and SW93-1015. (a) Volcano plots 
showing differentially abundant proteins. Blue dots at 24 and 72 hpi indicate proteins present at previous time points and grey dots indicate 
unique proteins for that time point, (b) Venn diagrams illustrating numbers of proteins with lower abundance (upper part) and higher 
abundance (lower part) in the apoplast of Desiree, Sarpo Mira and SW93-1015 from top to bottom at 6, 24, and 72 hpi from left to right. 



Secondly, we performed clustering by STEM, in order 
to find protein abundance patterns that occur during 
compatible and incompatible interactions in the apo- 
plast. Among 13 profiles with significant numbers of as- 
sociated proteins in Desiree, we found five profiles 
(number 26, 11, 9, 1, and 12 in order of significance) 
containing in total 156 proteins with lower abundance in 
the apoplast. In contrast, only one significant profile in 
Sarpo Mira (out of 10 profiles) and none in SW93-1015 
(out of 11 profiles) represented proteins with lower 
abundance (Figure 6a). In addition, one of the profiles 
representing proteins with higher abundance (49), was 
only significant for SW93-1015 and Sarpo Mira, and 
represented a large number of proteins for these clones, 
203 and 91 proteins respectively (Figure 6a). 

In order to find putative effector-targeted proteins that 
are differentially regulated in compatible and incompat- 
ible interactions, we selected candidates from profiles 
containing proteins with increased levels in the two in- 
compatible interactions (profile 39, 40, 42, 47, 48, and 49 
from Sarpo Mira, and profile 29, 39, 41, 42, 47, 48, and 
49 from SW93-1015) and compared them with the pro- 
teins associated with profiles representing decreased 
levels in the compatible interaction (profile 1, 9, 11, 12, 
26) (Figure 6a). We found 14 proteins from both incom- 
patible interaction profiles with increased abundance 
that showed decreased accumulation in the compatible 
interaction (Figure 6b; Additional file 5: Table S4). These 
included a lyase, a pectinesterase, a Kunitz trypsin in- 
hibitor and a dehydratase. In addition, we found 28 pro- 
teins with decreased abundance profiles in Desiree 
(Profile 1, 9, 11, 12, 26) but with increased abundance 



profiles in SW93-1015 (Figure 6c; Additional file 5: 
Table S4). Eight proteins with decreased abundance pro- 
files in Desiree were found in the increased abundance 
profiles of Sarpo Mira (profile 39, 40, 42, 47, 48, and 49; 
Figure 6d; Additional file 5: Table S4). 

Discussion 

We have carried out time course studies of one compat- 
ible and two incompatible interactions between P. infes- 
tans and potato by global transcriptomics and apoplastic 
proteomics using potato clones Desiree, Sarpo Mira, and 
SW93-1015. In response to P. infestans infection, there 
was a large overlap of transcripts and proteins with 
changed abundance irrespective of compatibility, which 
can be seen as a common signature for challenge with 
Phytophthora (Figures 2b and 3b). Among the differentially 
abundant proteins found in the apoplast, 50% were also 
differentially expressed at the transcript level (Figure lc). 
This is in line with previously reported data from yeast 
that showed that regulation of mRNA explains around 
40% of the protein concentration changes [39], and dem- 
onstrates the value of measuring both transcript and 
protein levels. 

We identified more than 1500 proteins in the apoplast; 
almost half of these are associated with 400 conserved 
functional domains (Figure la, Additional file 4: Table 
S3). The wide range of functional categories among apo- 
plastic secretome suggests a high complexity in this pro- 
tein fraction. Among the identified proteins 30% had a 
secretory signal predicted by TargetP, which is consistent 
with other apoplastic proteomics studies and suggests 




AN et al. BMC Genomics 2014, 15:497 
http://www.biomedcentral.com/1471-2164/15/497 



Page 9 of 18 



Table 2 Candidates for hypersensitivity initiation from microarray data 




RNA transcripts 




Increased in 


incompatible clones 


Gene ID 


Protein ID 


Name 


6h 


24 h 72 h 


PGSC0003DMG401 01 0883 


PGSC0003DMP400019223 


R2r3-myb transcription factor 


* 




PGSC0003DMG402004331 


PGSC0003DMP400007701 


Conserved gene of unknown function 






PGSC0003DMG400030212 


PGSC0003DMP400052593 


Nitrate reductase 


# 




PGSC0003DMG400001855 


PGSC0003 DMP40000332 1 


Beta-amylase PCT-BMYI 


* 




PGSC0003DMG400036101 


PGSC0003DMP400058205 


Gene of unknown function 


* 




PGSC0003DMG400019518 


PGSC0003DMP400033920 


Pseudo response regulator 


# 




PGSC0003DMG400011502 


PGSC0003DMP400020369 


PEP carboxylase kinase 


* 




PGSC0003DMG400020209 


PGSC0003DMP400035084 


Nodulin family protein 






PGSC0003 DMG40000 1 460 


PGSC0003DMP400002648 


Ninja-family protein Os03g0419100 






PGSC0003DMG400025240 


PGSC0003DMP400043812 


Axi 1 protein 


* 




PGSC0003DMG400020498 


PGSC0003DMP40003561 1 


lnositol-1,4,5-triphosphate-5-phosphatase 




* 


PGSC0003DMG400023702 


PGSC0003DMP400041014 


Calmodulin binding protein 




* 


PGSC0003DMG400023949 


PGSC0003DMP400041401 


Abscisic acid receptor PYL4 




# 


PGSC0003DMG400001771 


PGSC0003DMP400003164 


Dead box ATP-dependent RNA helicase 




# 


PGSC0003DMG400001387 


PGSC0003DMP400002507 


Ocs element-binding factor 




* 


PGSC0003DMG400010759 


PGSC0003DMP400019038 


Gene of unknown function 




* 


PGSC0003DMG400000731 


PGSC0003DMP400001424 


Response to dessication RD2 




* 


PGSC0003DMG400009982 


PGSC0003DMP400017630 


BTB/POZ domain-containing protein 




# 


PGSC0003DMG400003993 


PGSC0003DMP400007118 


Citrate binding protein 




# 


PGSC0003DMG400003057 


PGSC0003DMP400005490 


Osmotin 




* * 


PGSC0003DMG400015267 


PGSC0003DMP400026776 


Kunitz-type protease inhibitor 






PGSC0003DMG400026220 


PGSC0003DMP400045511 


Major pollen allergen Ory s 1 




* * 


PGSC0003 DMG400008 1 00 


PGSC0003DMP400014249 


KiTH-2 


# 


* * 


PGSC0003DMG400008099 


PGSC0003DMP400014248 


KiTH-2 




* * 


PGSC0003DMG400008098 


PGSC0003DMP400014247 


KiTH-2 


* 





Genes with increased expression in the two incompatible clones Sarpo Mira and SW93-1015, but not in compatible Desiree. Genes with a minimum 4-fold 
up-regulation in both Sarpo Mira and SW93-1015 at 6 or 24 hpi and stable or down-regulated in Desiree at 6 and 24 hpi were selected. Up-regulation is indicated 
with an asterisk. Bolded protein IDs and names indicate proteins for which abundance data was obtained. Symbol (*) indicates the time point in which the change 
was observed. 



that a large proportion of the secreted proteins are se- 
creted through unconventional secretory pathways [40]. 

At 6 hpi, a higher number of apoplast proteins increased 
in both incompatible interactions than in the compatible 
interaction, demonstrating a quicker response during the 
incompatible interaction (Table 1). This difference is even 
clearer at 24 hpi, when the number of proteins with in- 
creased abundance in the incompatible interactions is 
almost twice that of the compatible interaction. Another 
notable difference in the apoplast at this stage is that more 
proteins decreased in abundance in the compatible inter- 
action compared to the incompatible interactions. Inter- 
estingly, we did not observe this difference in the 
transcriptome, indicating a more pronounced suppression 
by effectors on secreted proteins than on transcription. At 
72 hpi, substantially more genes were induced in the 



compatible interaction than in the incompatible interac- 
tions. Most of these genes may have a role in metabolic 
and structural re-modeling due to infection. 

Although the total protein content is greater in the 
compatible interaction (Figure lb), fewer proteins show 
increased abundance. This suggests that the composition 
of the apoplastic secretome is important in the incom- 
patible interaction and that certain apoplastic proteins 
are suppressed by the pathogen during the compatible 
interaction (Figure lb; Tables 2 and 3). Major proteins 
that accumulated in the susceptible interaction were glu- 
canases and glucosidases (Additional file 2: Figure SI a) 
that might be involved in symptom development, as 
many more cells were visibly affected in the compatible 
interaction than in the incompatible interactions [34]. 
Evidence that the apoplastic secretome is altered during 
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Table 3 Candidates for hypersensitivity initiation from apoplast proteomics data 





Secreted proteins 


Increased in incompatible clones 


Protein ID 


Name 


6 h 24 h 72 h 


PGSC0003DMP400015021 


Cellulase containing protein 


* 


PGSC0003DMP400051976 


Glucan endo-1,3-beta-D-glucosidase 


* 


PGSC0003DMP400012140 


Polyamine oxidase 




PGSC0003DMP400056168 


Mitogen-activated protein kinase kinase kinase 




Q9LTJ3 


Putative uncharacterized protein At5g59350 




Q9SWI1 


Protein PH^OCHROME KINASE SUBSTRATE 1 


# 


DES_g15011_t1 


Uncharacterized protein 


# 


PGSC0003DMP400011041 


Leucine-rich repeat family protein 


# 


PGSC0003DMP400015631 


Aldo-keto reductase family 4 member G10 




PGSC0003DMP400027722 


Bacterial-induced peroxidase 




PGSC0003DMP400028029 


Alpha-glucosidase 




PGSC0003DMP400030201 


Strictosidine synthase 


* 


PGSC0003DMP400046981 


Kunitz trypsin inhibitor 


* 


PGSC0003DMP400049292 


Pattern formation protein 


* 


PGSC0003DMP400051822 


Actin-101 




PGSC0003DMP400001286 t 


Conserved gene of unknown function t 




PGSC0003DMP400002450 


Conserved gene of unknown function 




PGSC0003DMP400003176 


Zinc finger family protein 


* 


PGSC0003DMP400008097 


Rho guanine dissociation inhibitor 


* 


PGSC0003DMP400009086 t 


PR1 protein t 


# 


PGSC0003DMP400035459 


Triacylglycerol lipase 




PGSC0003DMP400040683 


ATP binding protein 




PGSC0003DMP400045856 


Pentatricopeptide repeat-containing protein 




PGSC0003DMP400057833 


Polyprotein protein 


* 


Q9FJT0 


Putative uncharacterized protein 


* 


Q9LDP1 


F28H19.2 protein (F2J6.15 protein) 


* 


Q9LVB8 


Probable carboxylesterase 120 (AtCXE20) (EC 3.1.1.1) 




Q9SCZ5 


Putative uncharacterized protein F2601 3.180 




PGSC0003DMP400005465 


Osmotin 


* * 


PGSC0003DMP400018074 


Phytophthora-inhibited protease 1 


* * 


PGSC0003DMP400021164 


Non-specific lipid-transfer protein 


* * 


DES_g20722_t1 


Uncharacterized protein 


* * 


022214 


Putative uncharacterized protein At2g41520 


* * 


064572 


Expressed protein (Uncharacterized protein) 


* * 


PGSC0003DMP400006538 


Aspartic proteinase nepenthesin-1 


* * 


PGSC0003DMP400012829 


Leucine-rich repeat receptor kinase 


* # 


PGSC0003DMP400026222 


Cytochrome P450 


* # 


PGSC0003DMP400038388 


Rhicadhesin receptor 


* * 


PGSC0003DMP400059150 


Gene of unknown function 


* * 


PGSC0003DMP400067598 


O-methyltransferase 


* * 


PGSC0003DMP400013560 


Gene of unknown function 


* * 


PGSC0003DMP40005391 1 


Beta-galactosidase 


* * 


PGSC0003DMP400012597 


Epidermis-specific secreted glycoprotein EP1 


* * 
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Table 3 Candidates for hypersensitivity initiation from apoplast proteomics data (Continued) 



PGSC0003DMP400016025 


Pentatricopeptide repeat-containing protein 








PGSC0003DMP400006168 


Glycosyl hydrolase family 


* 


* 


* 


PGSC0003DMP400017287 t 


Cytochrome P450 hydroxylase t 




* 


* 


PGSC0003DMP400035498 t 


Peroxidase t 






* 


PGSC0003DMP400064458 


Gag-pol protein 






* 


Q9FKV2 


Berberine bridge enzyme (FAD-binding and BBE domain-containing protein) 









Proteins with increased abundance in the secretome for the two incompatibie clones Sarpo Mira and SW93-1015, but not for compatible Desiree. Proteins with 
a minimum 4-fold increase in both Sarpo Mira and SW93-1015 at 6 or 24 hpi and stable or decreased abundance in Desiree at 6 and 24 hpi were selected. 
Up-regulation is indicated with an asterisk. The symbol t indicates proteins that were not identified using at least one peptide specific for the listed gene. Symbol 
(*) indicates the time point in which the change was observed. 



pathogen attack has previously been found in cell cul- 
tures [41], but such a large-scale perturbation of the host 
apoplastic secretome in planta is unprecedented. 

Little is known about the genes involved in develop- 
ment of HR in potato. There were few genes that were 
uniquely induced in both resistant potato clones; notably 
glutaredoxins, MYB transcription factors and a zinc fin- 
ger protein (Tables 2 and 3; Additional file 4: Table S3). 
Plants with impaired glutaredoxin activity are resistant 
to the necrotroph Botrytis [42] and SA-inducible glutar- 
edoxin is involved in suppression of the jasmonic acid- 
responsive genes in Arabidopsis [43], suggesting a role 
in hormone crosstalk during the defense response. Plant 
U-box proteins, also unique to incompatible interac- 
tions, is a gene family with diverse functions linked to 
salicylic acid and PAD4 mediated pathways [44]. WRKY 
transcription factors did not show this pattern and many 
were initially down-regulated, although WRKY8 has 
been implicated in defense-related MAPK signaling and 
P. infestans resistance in Nicotiana benthamiana [19]. 

Sarpo Mira, which carries five different resistance 
genes [45], showed relatively few changes after challenge 
with P. infestans. This might be due to specific R gene 



recognition and responses to P. infestans. Although the 
difference between specific responses in the two resist- 
ant genotypes was large, both still display full resistance 
against the pathogen and initiate similar macroscopic 
HR in detached leaf assays [34]. Our analysis suggests 
the possibility that pathogen recognition occurs differ- 
ently in Sarpo Mira and SW93-1015. Sarpo Mira might 
mainly depend on effector recognition by resistance pro- 
teins, while SW93-1015 may have additional recognition 
and response mechanisms. At 6 hpi in SW93-1015 two 
LRR receptor-like kinases are induced, perhaps leading 
to defense responses unique to SW93-1015. 

The fact that certain proteins in the apoplast might be 
selectively suppressed during the compatible interaction 
but not the incompatible interactions allows us to screen 
for putative effector targets. In this unique type of 
screening described here, we found several candidates 
that were single proteins from large gene families. In 
addition to the well characterized cysteine protease 14 
(PGSC0003DMP400017510), we found another cysteine 
protease (PGSC0003DMP400045977) as well as a cysta- 
tin homologue (a putative cysteine protease inhibitor 
PGSC0003DMP400018076). Two other proteins of 
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Figure 5 Gene expression and protein abundance in the apoplast of known effector target cysteine protease CI 4. Transcript and protein 
levels of the effector target cysteine protease C14 decreased at 6 hpi in the Desiree apoplast and at the same time was transitionally 
up regulated. 
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Table 4 Putative effector targets 



Protein ID 


Name 


Desiree 
6 h 


Desiree 
24 h 


PGSC0003DMP400017664 


Transcription factor 


# 




Q9SZ87 


Reverse transcriptase (RNA-dependent DNA polymerase) 


* 




Q9SPM5 


GDA1/CD39 (nucleoside phosphatase) family 


* 




Q8LBK6 


Glutaredoxin 


* 




Q9SSD1 


Leucine rich repeat 


# 




PGSC0003DMP400018076 


Cysteine protease 


* 




PGSC0003DMP400045977 


Cystatin 


* 




PGSC0003DMP400024264 


Conserved gene of unknown function 


* 




PGSC0003DMP400020961 


5'-nucleotidase sure 


* 




081459 


ATPase family associated with various cellular activities (AAA) 


# 


# 


PGSC0003DMP400049952 


Pectinesterase-2 


* 


# 


F4J1D9 


MutS domain II 




# 


PGSC0003DMP400008301 


Pto-interacting protein 1 


* 


# 


Q8VZM7 


Protein of unknown function (DUF1012) 


* 


* 


Q9FJV2 


FBD, Leucine Rich Repeat 


* 


# 


PGSC0003DMP400006604 


Beta-galactosidase t 


# 


# 


PGSC0003DMP400021388 


HSP transcription factor t 


* 


# 


PGSC0003DMP400001052 


Peptidyl-prolyl cis-trans isomerase 


* 


# 


PGSC0003DMP400054971 


Phospholipase A1 t 


* 


* 


PGSC0003DMP400065569 


Gene of unknown function 


* 


# 


DES.g31837.t1 




* 


# 


PGSC0003DMP400056405 


Histidine kinase 3B 




# 


PGSC0003DMP400016823 


Kunitz-type proteinase inhibitor 






PGSC0003 DMP400044843 


Serine-threonine protein kinase, plant-type 


* 


* 


PGSC0003DMP400004370 


Glutamate decarboxylase 4a 


* 


# 


Q9FGS4 


Fe-S metabolism associated domain, Quinolinate synthetase A protein 


* 


* 


A0ME24 


Putative uncharacterized protein 


* 


# 


081861 


Glycosyl hydrolases family 18 


* 


# 


rojL-UUUjUIVIr4UUUUo/ ID 


Subtilase family protein 




# 


Df^c f'nnn? r~\k ad a r\nr\£.Q i o c 
rbbLUUU3UMr4UU(JoojoD 


Conserved gene of unknown function 




* 


Drcmnn^ ni\ ad a r\r\nn/\ <;n 


Beta-galactosidase 


* 


* 


rUbLUUU3UIVIr4UUUDD43y 


Protein kinase atmrkl 






MUbLUUU:5UMr4UUU I / j I U 


Cysteine protease 14 


* 


* 


p(^rnnn3ni\/iP4onrn 771 q 

r jluuu j l^i vir ilaju i / / i ? 


Ml tirnnnpr rwiHacp 




# 


Q8VZG2 


Domain associated at C-terminal with AAA, ATPase family associated with various cellular activ- 
ities (AAA) 




# 


PGSC0003DMP400048051 


Glutathione S-transferase T5 




# 


Q4FE22 


Potato inhibitor 1 family 




# 


PGSC0003DMP400049086 


Thioredoxin peroxidase 




# 


Q9FMR7 


Mitochondrial fission ELM1 




# 


PGSC0003DMP400025698 


High mobility group protein 




# 



Proteins with decreased abundance in the secretome for the compatible clone Desiree, but not for the two incompatible clones Sarpo Mira and SW93-1015. 
Proteins with a minimum 8-fold down-regulation in Desiree and stable in both Sarpo Mira and SW93-1015 at 6 or 24 hpi were selected. 8-fold down-regulation is 
indicated with an asterisk. The symbol t indicates proteins that were not identified using at least one peptide specific for the listed gene. Symbol (*) indicates the 
time point in which the change was observed. 
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(See figure on previous page.) 

Figure 6 Identification of candidates for effector targets, (a) Significant profiles from STEM clustering of apoplastic proteins at 0, 6, 24, and 72 

hpi in the three potato clones. Protein intensities from four replicates for each sample were used for clustering. Only proteins with statistically 
significant differences for at least at one time point were used. Each box represents a protein abundance pattern during the time course. Top left 
of each box is the profile number defined in STEM and bottom left of each box indicates the number of proteins that fit the defined profile 
pattern, (b) Proteins found with decreased abundance profiles in Desiree, but increased abundance profiles in both Sarpo Mira and SW93-1015. 
(c) Proteins found with decreased abundance profiles in Desiree but with increased abundance profiles in SW93-1015. (d) Proteins found with 
decreased abundance profiles in Desiree but with increased abundance profiles in Sarpo Mira. 



interest were a Kunitz-type proteinase inhibitor (PGSC 
0003DMP400016823) and a subtilase family protein 
(PGSC0003DMP400008705) that also was identified in 
this screen. With the exact identity of these genes it is 
now possible to investigate these functionally. 

Conclusion 

We carried out a time course study of potato transcrip- 
tomic and proteomic responses after P. infestans inocu- 
lation, both in compatible and incompatible interactions. 
Studying samples from the same time course for both 
gene expression and protein abundance lead to identifi- 
cation of potential targets for plant proteins directly ma- 
nipulated by the pathogen. By using two phenotypically 
different resistant potato clones we identify several tran- 
scripts and proteins that only show increased expression 
or abundance in both incompatible potato clones that 
are potentially involved in resistance. In summary, we 
provide ample number of transcripts and proteins that 
could be used in targeted studies on one of the most 
agriculturally important plant-pathogen interactions. 

Methods 

Plant material and growth conditions 

Three potato clones Desiree, Sarpo Mira and SW93- 
1015 were used. Sarpo Mira and SW93-1015 are two 
highly resistant potato clones with slightly different re- 
sistance reactions [34]. Sarpo Mira is a potato cultivar 
that recognizes five different effectors from P. infestans 
[45] and shows a classical HR reaction in response to P. 
infestans inoculation. SW93-1015 is a breeding potato 
clone which is consistently resistant to Swedish P. infes- 
tans populations with reduced HR expansion and a weak 
cpr genotype [34]. Plants were grown in a growth cham- 
ber with controlled conditions set at 20°C with a 16:8 
lightdark cycle and 70% relative humidity. Five-week-old 
plants were transferred to an infection chamber with 
100% humidity and 10:14 light:dark cycle. After 6 hours, 
plants were sprayed with an encysted zoospore suspen- 
sion from P. infestans isolate SE-03058 [46] until the leaf 
surfaces were fully saturated with the zoospore suspen- 
sion (15000 sporangia/ml). Relative humidity was main- 
tained at 100% for 2 days after inoculation and then 
adjusted to 90% for the rest of the experiment. Samples 
were collected at 6, 24, and 72 hpi. For RNA and 



apoplastic secretome sampling, fully-expanded upper 
leaves were collected. 

RNA isolation and microarray analysis 

Samples were collected from three independent bio- 
logical experiments. Four to six leaves frozen in liquid 
nitrogen were homogenized to a fine powder using 
FastPrep-24 (MP-bio, Santa Ana, USA) with 2 mm 
beads. Leaf tissue (100 mg) was weighed out and put 
into RNase-free tubes before extraction. RNA extrac- 
tions were performed using the RNeasy Plant Mini kit 
(Qiagen GmbH, Hilden, Germany). Samples were DNase 
treated and cleaned using the Qiagen RNA cleanup kit. 
RNA concentration and purity (260/280 nm > 1.8) was 
checked by a ND-1000 NanoDrop (Wilmington, USA) 
and integrity of the samples were analyzed with an 
Experion™ Automated Electrophoresis System (Bio-Rad 
Laboratories, Hercules, USA). For mRNA expression 
analysis, a custom-made Agilent expression array (JHI 
Solanum tuberosum 60 k vl; ArrayExpress accession A- 
MEXP-2272) based on predicted transcripts of the Sola- 
num tuberosum genome (version 3.4) was used according 
to the supplier's directions (One-Color Microarray-Based 
Gene Expression Analysis Low Input Quick Amp Labeling 
v. 6.5; Agilent). Data were extracted from each scanned 
array image using Feature Extraction software (v. 10.7.3; 
Agilent). The array data is deposited in ArrayExpress: 
E-MTAB-1515. 

Apoplastic secretome sample preparation 

Three biological replicates of un-infected control sam- 
ples were collected. For each apoplast sample, four fully 
expanded middle leaves from two plants were subjected 
to independent apoplast isolations. In order to collect 
un-infected control samples, plants were transferred to 
the humid chamber and kept in the chamber for 6 hours 
under the same experimental conditions as for infection 
experiments. Apoplast isolation was performed by using 
vacuum infiltration as described [34]. Apoplastic fluids 
collected from two plants for each sample were pooled, 
dissolved in 6x SDS-PAGE buffer containing DTT, and 
denatured at 95°C for 3 minutes. Pooled samples (30 ul 
of each) were loaded on polyacrylamide gels and sepa- 
rated for 2 cm with SDS-PAGE. After staining with 
Coomassie, the gel lane from each sample was cut into 
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about 1 mm 2 pieces. Each lane was kept separate from 
this point on, generating five subreplicates for each con- 
trol sample and four subreplicates for each infected sam- 
ple. Samples were then subjected to in-gel tryptic 
digestion with incubation (modified sequencing grade; 
Promega, Madison, WI, USA) overnight at 37°C. Pep- 
tides were extracted in 50-80% acetonitrile and excess 
acetonitrile was vaporized using centrifugation under 
vacuum. De-salting was performed using UltraMicro 
spin columns (Nest group). The whole experiment was 
repeated twice. 

Mass spectrometry 

MS analysis was performed on a LTQ Orbitrap XL with 
an Eksigent nano-LC system (Eksigent technologies, 
Dublin, CA, USA). A 5 ul sample was injected and sepa- 
rated at a flow rate of 300 nl/min with a 90 minute 
gradient. The four most intense ions were selected in 
data-dependent mode for fragmentation in the linear ion 
trap, for details see [47]. Files were converted to mzML 
[48] and Mascot Generic Format (MGF) using Proteo- 
Wizard [49] and uploaded to the Proteios Software En- 
vironment [50]. MGF files were used for MS/MS 
identification, and mzML files for feature detection. 
Identification searches were performed in Mascot (www. 
matrixscience.com) and XITandem (www.thegpm.org/ 
tandem) in a database consisting of all Solanum proteins 
in UniProt (http://www.uniprot.org), all annotated pro- 
teins from the potato genome project (http://www. 
potatogenome.net, [6]) and predicted proteins from de 
novo assembled transcripts from three potato clones De- 
siree, Sarpo Mira and SW93-1015, plus an equally sized 
decoy part consisting of the reversed protein sequences 
(in total 670584 proteins). The MS mass tolerance was 
set to 5 ppm and MS/MS fragment tolerance to 0.5 Da, 
with one missed cleavage allowed. Cysteine carbamido- 
methylation was set as fixed and methionine oxidation 
as variable modification. FDR was subsequently calcu- 
lated for the combined search results at the peptide- 
spectrum level, and filtered at a FDR of 1% as described 
previously [50,51]. To quantify possible peptides, msln- 
spect [52] feature detection was performed on the 
mzML files from Proteios using default settings. The fea- 
tures were matched to identifications with a retention 
time tolerance of 0.2 minutes and an m/z tolerance of 
0.005 Da as well as a requirement of same charge and 
LC-MS/MS run. To facilitate the differential expression 
analysis and to propagate sequences between the runs, 
alignment was performed using our recently developed al- 
gorithm within Proteios [47,53]. A report of the peptide 
features corresponding between runs was exported for 
further analysis. Peptide data used for further analysis is 
found in Additional file 6: Table S5 and LC-MS runs 
are listed in supporting information, Additional file 7: 



Table S6. The mass spectrometry proteomics data have 
been deposited to the ProteomeXchange Consortium 
(http://proteomecentral.proteomexchange.org) via the PRIDE 
partner repository with the dataset identifier PXD000435 
and DOI 10.6019/PXD000435. 

Bioinformatics and statistical analysis 
RNA-sequencing 

We obtained approximately 52 million reads for each 
potato clone from RNA samples collected at 24 hpi 
using paired-end libraries from Illumina HiSeq 2000 ma- 
chines. Clean paired-end reads were then de novo as- 
sembled using Trinity (version r201 1-1 1-26) to build 
transcript contigs. Protein coding sequences from these 
contigs were extracted by gene prediction program Au- 
gustus (version 2.5.5) that was used for constructing a 
proteomics database. BLASTX was used to obtain anno- 
tations from the potato genome. The RNA-seq data has 
been deposited in ArrayExpress: E-MTAB-1712. 

Microarray 

Probe intensities were background-corrected and nor- 
malized using the quantile method in the Limma R 
package [54]. Genes with p-values below 0.05 after ad- 
justment with the Benjamini-Hochberg method were 
regarded as significant. The projected GO annotation 
(Additional file 8: Table S7) [55] for the differentially 
expressed probes was then analyzed for enrichment of 
Gene Ontology (GO) terms using GOEast [56] with de- 
fault settings (adjusted p-value < 0.1). 

Apoplastic secretome 

Peptides with a FDR of < 0.01 were selected for further 
analysis. For normalization, we used the Eigen MS 
method incorporated in DanteR (v0.2) that uses Eigen 
values to find trends in the data for normalization. Only 
peptides that were identified in at least two of the subre- 
plicates of any sample were selected and normalized. In 
DanteR, data was filtered, missing values imputed and 
protein level intensities for the leading proteins calcu- 
lated using Q-Rollup function with 2 to 3 peptides per 
protein based on median peptide intensities [57-59]. 
Proteins with single peptide hits were included for fur- 
ther analysis. Differentially expressed proteins were cal- 
culated using one-way ANOVA for all time points from 
each potato clone with their relative controls. Fold 
change estimates were performed based on linear model 
comparison of each time point with the relative control 
for each potato clone. After Benjamini-Hochberg (BH) 
p-value adjustment, differentially expressed proteins 
with adjusted p-value < 0.05 were selected for further 
analysis. Default settings in STEM [60] were used and. 
Standalone InterProScan [35] with default settings was 
used to find conserved domains and peptide fingerprints 
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for apoplastic secretome protein classification [61]. Tar- 
getP was used for secretory signal peptide identification 
[36]. MapMan analysis was performed to explore gene 
pathways based on sequence identity to Arabidopsis 
genes [61]. 

Additional files 



Additional file 1: Table SI. Clone-specific peptide identifications. ORFs 
were predicted based on RNA-seq data in Desiree, Sarpo Mira and 
SW93-1015. 

Additional file 2: Figure SI. Total peptide intensity for proteases and 
peptidases, a, Sum of total intensities of 49 proteins in the apoplast 
identified as glucanases and glucosidases (4 replicates) b, Sum of total 
intensities of 78 proteins in the apoplast identified as proteases and 
carboxypeptidases (4 replicates). Error bars indicated ± SD of the mean. 
Figure S2. Functional analysis of gene expression by MapMan. Bins for 
genes related to hormones (a) and transcription (b) for Desiree, Sarpo 
Mira and SW-1015 at 6, 24 or 72 hpi with an adjusted p-value of <0.05 
are shown. 

Additional file 3: Table S2. List of enriched Gene Ontology terms 
across the time course. Enrichment was done with GOEast using default 
settings and a cut-off of p < 0.1 . 

Additional file 4: Table S3. Gene expression during Phytophthora 
infection. Microarray data was normalized and analyzed with the Limma 
R-package. Expression ratios and corresponding adjusted p-values are 
given for all microarray probes with matching potato transcripts. 

Additional file 5: Table S4. Protein identifiers for resistance candidates 
identified in STEM profiles. Excel file with identifiers correspond to 
profiles presented in Figure 6b-d. 

Additional file 6: Table S5. Peptide data during Phytophthora infection. 
Raw mass spectrometer peptide intensities. 

Additional file 7: Table S6. LC-MS file names. Raw data was have been 
deposited to the ProteomeXchange Consortium (http://proteomecentral. 
proteomexchange.org) via the PRIDE partner repository with the dataset 
identifier PXD000435 and DOI 1 0.601 9/PXD000435. 

Additional file 8: Table S7. Gene ontology terms for the potato 
genome and custom-made Agilent expression array. The file (Extended_ 
GOterms_JHI_Solanum_tuberosum_60kv1_expression_array.txt) is 
formatted for use in GOEast [57]. 
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